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We show that radiative coupling between two multilevel atoms having near-degenerate states can 
produce new interference effects in spontaneous emission. We explicitly demonstrate this possibility 
by considering two identical V systems each having a pair of transition dipole matrix elements 
which are orthogonal to each other. We discuss in detail the origin of the new interference terms 
and their consequences. Such terms lead to the evolution of certain coherences and excitations 
which would not occur otherwise. The special choice of the orientation of the transition dipole 
matrix elements enables us to illustrate the significance of vacuum induced coherence in multi-atom 
multilevel systems. These coherences can be significant in energy transfer studies. 
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A multilevel atom having closely lying energy states (energy separation of the order of the natural linewidth) can 
show interferences in the decays from those closely lying levels to a common ground state. This is due to the fact 
that, both the decay channels are coupled via the same continuum of the vacuum, creating the interfering path ways. 
The resulting coherence in the system is known as vacuum induced coherence (VIC). Occurrence of this coherence 
requires a stringent but achievable condition - i.e. the transition dipole matrix elements involving the decay processes 
should be non- orthogonal The manifestation of VIC in atomic systems has given rise to a myriad of fascinating 
phenomena [1-17]. All these studies deal with a single multilevel atom or equivalently with an ensemble of non- 
interacting multilevel atoms (e.g. very low density atomic gas systems). However, VIC in coupled atomic systems has 
remained unexplored. In this paper, we consider the role of VIC in two radiatively coupled multilevel atoms. 

We start by recalling some of the consequences of VIC in a single atom. It was first shown by Agarwal W that 
""^ : population gets trapped in degenerate excited states of an atom due to interference in decays channels. Recently, 
there has been renewed interest in this subject particularly in the context of coherently driven systems [§-^|. Harris 
and Imamoglu were the first to discover the possibility of achieving lasing without population inversion in systems 
where two excited states were coupled to a common continuum |(| (See also @,&). It has also been been observed that 
narrowing of spontaneous emission can be obtained by making use of the VIC |3|,|9| • Quantum beat has been observed 
in spontaneous emission which showed pronounced beat structures determined by the energy separation of the closely 
lying states |l(],[Il|]. The VIC also leads to cancellation of spontaneous emission [ p^ - |l4"||. Zhu and coworkers have 
experimentally demonstrated the quenching of spontaneous emission in sodium dimers |12|. Further Scully, Zhu and 
coworkers proposed many schemes with different configurations demonstrating the possibility of obtaining quenching 
of spontaneous emission It was also reported that in the presence of VIC, the resonance fluorescenc e |15|]16[ 

and other spectral line shapes || become sensitive to the phase of the control laser. Knight and coworkers |lf|] have 
demonstrated the possibility of controlling spontaneous emission by varying the relative phase of two control lasers in a 
four-level system. The question of requirement of non-orthogonal dipole moments has been addressed and alternative 
possibilities have been suggested [ pd] , p^|Jl^Jr7[ ] . 

All the above works [2-17] correspond to a single atom system, or equivalently an ensemble of non- interacting atoms 
where the average inter-atomic distance is much larger compared to the wavelength of the emitted radiation. However, 
when the inter-atomic distance becomes comparable to the wavelength, the dipole-dipole (dd) coupling between the 
atoms via vacuum gives rise to collective effects. Our usage of dd interaction should be understood in the sense of 
retarded (and complex) dipole-dipole interaction. The classic example is Dicke superradiance |T§| ], where the atoms 
in their excited state decay much faster compared to that of the single atom case. The collective effects in atoms have 
been extensively studied [18-30]. Recently experiments have been reported to observe collective behavior with two 
identical atoms |l!],|2(| . The dd interaction has been shown to produce two-photon resonance |^l[ and frequency shifts 
in emission |22] . The energy exchange between two coupled systems is discussed in |23| . Many interesting features of 
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dd interaction in the context of atoms interacting with a squeezed vacuum EJ] , and inside bandgap materials |25| have 
been reported. Mayer and Yeoman 26 have considered two-atom laser in the presence of the atom-atom interaction. 
Quantum jump from two dipole interacting V-systems giving rise to new fluorescence periods has been reported by 
Hegerfeldt and coworkers P7|| . Meystre and coworker s pq| found that the dd interaction leads to the occurrence of 
dark states in the fluorescence of two moving atoms |29|]. Finally note that the local field effects in a dense media 
are also a consequence of dd interaction J30|. All the dd interaction related effects can be understood as due to the 
exchange of virtual photons between the atoms. Most of the existing literature concerns two-level atoms. 

In this paper, we consider two identical V-systems having two closely lying excited states (as shown in Fig. 1). 
The two atoms get coupled by the exchange of radiation. We would specifically show how the radiative coupling 
in multilevel systems can lead to a population transfer from \1a) to \2b) even if the corresponding dipole matrix 
elements are orthogonal. 




Atom A 




Atom B 



R 

FIG. 1. The two identical V systems under consideration. The distance between the two atoms is R. The transition dipole 
matrix elements di and d^ are chosen to be orthogonal with each other. The energy separation between the excited states is 
T18. We discuss the situation where initially atom A is in excited state |1a) and atom B is in ground state |3s). Possibility of 
dd interaction induced excitation in the |3s) — ► |2s) transition is indicated. 



The organization of the paper is the following: In Sec. II we derive the equations for the dynamical evolution of 
the two V-systems in master equation formalism. In Sec. Ill we interpret the terms appearing in the master equation 
and discuss the dependence of these terms on the geometry of the atoms in detail. In Sec. IV we present the numerical 
results on the evolution of the excited state coherence and the excitation probabilities. In Sec. V we present results 
for the case of magnetic degeneracies and we also discuss how the new coherence effect can be monitored. In Sec. VI 
we present the concluding remarks. 



II. DYNAMICAL EVOLUTION OF TWO V-SYSTEMS 



We consider [Fig. I] two identical V-systems (say A and B) in free space, having two near-degenerate excited states 
1 1 M ) and \2fj) (fj, = A, B) with the level separation %8. The ground states of the atoms are represented by |3 M ). Let 
w\ and u>2 be the atomic frequencies corresponding to 1 1 M ) <-> 1 3 M ) and |2 M ) <-> |3 P ) transitions respectively. Let the 
position vectors of the atoms be xa and xb- Both the atoms couple with the vacuum field, which is given by 

£„(f) = £( + '(f) + £(-)(£); (1) 

where E^ + ' (E^~') represents the positive (negative) frequency part the vacuum field at x and is defined as 

= ^(^) 1/2 a ks i ks e^, (2) 

k,s 

and E(~} is the Hermitian conjugate of E^ + \ Here aks (^L) represent the annihilation (creation) operator of a field 
mode having propagation vector k and polarization s; uJk(= kc) represents the angular frequency corresponding to 
the fc-mode. The corresponding unit polarization vector is denoted by tks- 

In what follows we consider the following physical process: Initially atom A is taken to be in excited state 11a) and 
the second atom in ground state \2>b)- To highlight the effect of new interference terms, we specifically consider the 
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case when the two transition dipole matrix elements d\ and d 2 are orthogonal to each other. In the absence of atom 
B, the VIC cannot be created in the excited states of atom A because the transition dipole matrix elements d\ and 
o?2 are orthogonal. However, the radiative coupling can lead to evolution of the excited state coherences. We will also 
show a manifestation of this coherence in the dynamical evolution of the atomic population. 
The total Hamiltonian for the atoms and the field system is given by 

H = H a + H f + H I , (3) 

where the unperturbed Hamiltonian for the atoms and field are 

H a = h + ^2/3^) , (4) 

li=A,B 



Hf = ^ huksd{ s a ks 



(5) 



and the interaction Hamiltonian is 



H, 
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(6) 



■ Y, [djf ] ■ E v (x^) + /i.e.] 0* = A, B). 
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The atomic transition operators introduced in (Q), are given by 



it =: 



(7) 



Thus dp, ajt fijj represent the atomic lowering and raising operators respectively corresponding to the |3 P 
(|2p) <-> |3 M ) ) transition. The dipole matrix elements are represented by 



(8) 



For simplicity we consider a situation where the transition dipole matrix elements of atom A are parallel to the 
transition dipole matrix elements of atom B 



4 A) and d 2 r \\d 23 



(9) 



We also assume that 



4 A) -4 A) * = o. (io) 

Thus the index /i in the right hand side of Eq. (@) can be dropped and the dipole matrix elements can be rewritten as 



dV=d 3 \^){^\. 



(11) 



Here we note that dj, in general, can be complex (see for example Sec.V). Using Eq. (|ll|), the interaction Hamiltonian 
in (|6|) reduces to 



Hi = - Yj [{di&l + ■ E v {x^) + h.c. 



V=A,B 



(12) 



We work in the interaction picture by transforming ( Jl2] ) 

Hi(t) = e ^+ H ^Hie-^+ H ^ 

Y [(didte^ + ^te^*) • (iW(i (11 () + i(-'(2 il) i)) + 



ti=A,B 



(13) 
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where, 
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1/2 



Ok. 



Ase^"-^ and E^fo,*) = [£<+>(£ M ,t)]t. 



(14) 



Let the density operator of the combined atom-field system in the interaction picture be represented by p(t) which 
satisfies the Liouvillc equation of motion 



dp 
dt 



[Hi(t),p] 



(15) 



To derive useful information about the evolution of the atomic system, we derive a master equation for the reduced 
atomic operator by using the standard projection operator techniques [Q. We make certain approximations in deriving 
the master equation: (a) at t = 0, p(0) can be factorized into a product of atom [p a ] and field [pf] density operators, 
i.e. p(0) = p a (0)pf(0). Furthermore, we invoke (b) the Born approximation and (c) the Markoff approximation. The 
Born approximation depends on the weak coupling between the vacuum and the atoms. The Markoff approximation 
holds because the vacuum has fairly flat density of states. Using the above approximations and tracing over the field 
states, the density matrix equation for the atoms becomes 



dp a 
dt 



-2 lim / drTrfiH^lH^t-r^pf^pa}] 



(16) 



The trace over the field operators inside the integral in Eq. ( jlif ) is calculated using the following relations: 

Tr /(P/ a L«fc's') = °> Tr /(P/ a fc s aL') = Skk'S ss ', 

Tr/ (p f a ks ak's') = Tr/(p/4 s 4, s ,) = 0. (17) 

One also uses the rotating wave approximation (RWA) to drop the anti-resonant terms like a^cij,, /3^i/3 M and 

PuPu m (0)- Using the above conditions and carrying out a long algebra, we obtain the master equation for the 
atomic density operator 



dp 

dt 
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where, 
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Here, R = S A — x B . Since the states |1 M ) and |2 M ) are closely lying, we have set uj\ = u) 2 = u>o- The sufSx a of p a in 
Eq. (|l^) has been dropped for brevity. We have also dropped the Lamb shift terms associated with the spontaneous 
emission of the individual atoms. The summation over the polarization components is evaluated using the relation 



^^(efcs)^t (e^ )^ — Spy k^kv, 



(20) 



where ki represents the direction cosine of k/\k\ along the I th Cartesian axis. Taking the limit V — > 00 and replacing 
the summation over k by integration over the continuum of the field modes, the terms in (|l9| ) become 



7i 



_ 2\dj\ 2 /woy 
3h lei 



T. t = i {dt • Im X -d*), ^ = \ (di ■ Re X -d*), 
T vc = - (d 2 ■ Im X ■(%), il vc = - (d 2 ■ Re X 



(21) 



where X is a tensor whose components are given by 



tJ 2 



dx Au dx Al 



R 



(k = cjo/c) 
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e lfcoi? ;i?= \x A -x B \ 



(22) 



III. INTERPRETATION OF DIFFERENT TERMS IN THE MASTER EQUATION (18) 

The 27^ in Eq. ( f2l| ) represents the single atom spontaneous decay rate from the state \i) to the state |3). Rest of the 
coefficients in ( |l8| ) are related to the coupling between the two V-systems. This coupling is produced by the exchange 

of a photon between the two systems. The dipole-dipole interaction manifests itself through the tensor X defined by 

Eq. (p2|). The tensor X a p {%A, %B>u) has the following meaning: It represents the a th component of the electric field 
at the point x A , produced by an oscillating dipole of unit strength along the direction and located at the point 
xb pl[ . In the limit c — * 00 (fco — ► 0), it reduces to the static dipole-dipole interaction. I\ and Oj represent the dd 
couplings which are related to the decay and level shift of the collective atomic states. These coefficients couple a 
pair of parallel dipoles and are well known [Q], particularly in the context of collective effects in two level atoms. The 
new coherence terms T vc and il vc are the dipole-dipole cross coupling coefficients, which couple a pair of orthogonal 
dipoles. The meaning of such terms can be clearly understood by calculating the evolution of the population in the 
state \3 A , 2b), given the initial condition \l Al 3b)- From the master equation (Eq), one can show that at t ~ 0, 



and hence 
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^(U, 3 B \ P \3 A , 2 B ) = -(Kc + M* vc )e m , (23) 



(3 A , 2 B \p\3 A , 2 b ) = -(T„ c - iQ vc )e- i5t {l A , 3 B \p\3 A , 2 B ) 
at 



-{T* vc +iQl c )e^{3 A ,2 B \p\l A ,3 B ) 

-2\T vc - l n vc \ 2 ( S ^). (24) 



Thus to the lowest order in T vc and £l vc , we obtain 



,2 

\2 



(3 A , 2 B \p\3 A , 2 B ) « 4|r wc + iSl vc f ^ sm . (25 ) 

Therefore the radiative process, in which atom A in the excited state \1 A ) loses its excitation which in turn excites 
atom B to the state |2b), is possible only because of T vc and Q, vc terms in the master equation (O). Note further 
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that such terms start becoming insignificant as US - the energy separation between the two excited states increases. 
Such interference terms occur in the master equation even when the transition dipolc matrix elements d\ and d 2 
are orthogonal. Such contributions come from the second term in Eq. (|22]). In the rest of the paper, we study in 
detail various consequences of these interference terms which could be large and could significantly contribute to the 
dynamics of the system when the atomic separation is smaller than A. 

Let us consider a geometry where atom A is placed at the origin of a Cartesian co-ordinate system and the position 
vector of atom B is R (as shown in Fig. 2). R makes an angle 9 with the z axis. Let us assume d\ = xd and di = yd. 
All the radiative coupling terms in the master equation can be written down explicitly as 



(26) 
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y(Pi - sin 2 9 cos 2 <j) 
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z ^ Re Xyx = 


37 n 

= sin 9 sin 6 cos ( 
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1> Qr] 



where <f> is defined as in the Fig. 2, and 



R, - - , ■ Q, - 3 " - £ 2 3 ^3 ' 

cos C sin C, 



cos ( 


sin C 


cosC 




c 


c 2 


c 3 ' 


Q,= c ° c sC 


sin C 


cosC 


sinC 






h «• 


c 3 ' 





^ . ^ , C = fcoi?. (27) 



Z A 




FIG. 2. The geometry under consideration where the dipole matrix elements di and are taken to be real and aligned 
along the x and y directions respectively. 

In the following, we examine the behavior of the cross coupling coefficients responsible for the new coherence effects 
in different geometries. In Fig. 3 we plot these coefficients as a function of the distance between the two atoms. In 
Fig. 3(a) we have plotted Tx and r„ c , and in Fig. 3(b) we plot Six and £l vc , for comparison. Clearly, the values of the 
cross coupling coefficients are comparable with the Ti and Qi values. The value of £l vc becomes significantly large for 
R < A/2. However for R — * 0, the terms £li,tt vc diverge, whereas the terms r.i,r„ c — > 1. 
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Further, in Fig. 4 we examine the atomic position dependences of these coefficients. We have plotted the coupling 
coefficients as a function of </>. Here we have fixed 9 = tt/2; i.e. both the atoms are lying in the xy-plane. Again for a 
comparison, we have plotted Ti and r„ c in Fig. 4(a), and Sli and Cl vc in Fig. 4(b). We observe the following special 
cases: 

Case I: If 9 = nn, then T vc = Q vc = 0; i.e. if R is perpendicular to the plane containing d\ and d,2, the interference 
terms in the master equation drop out. 

Case II: When (f> = mr/2, the coherence terms T vc = Q vc = 0; i.e. when the second atom is placed in a position such 
that R is along one of the dipoles d\ or e?2, then again the interference terms drop out. Thus the interference effects 
in the radiativcly coupled systems are sensitive to the geometry. 
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FIG. 4. The dd-coupling terms as a function of the azimuthal angle <j)- Here 6 = n/2 and the atomic separation is taken to 
be A/4. All coefficients are scaled with 7. 

IV. NUMERICAL RESULTS 

In this section we present the numerical results that demonstrate the effect of the interference terms on the dynamics 
of the radiatively coupled multilevel systems. We use fifth order Runge-Kutta method for the numerical solution of 
the master equation ([l8|). For numerical solutions we use the initial condition that at t — 0, the first atom is in excited 
state 1 1^4) and the second atom is in the ground state \3b)- 

In Fig. 5, we have plotted the density matrix element p$ = 0-A\(^B\p(t)\2A)\^B) which represents the coherence 
in the excited states of atom A when atom B is in ground state |3s). It is clear from Fig. 5 that the interference 
terms r„ c , Q. vc in the master equation result in finite coherence in atom A. Otherwise, when T vc — f2„ c = 0, such 
coherences vanish. It is important to note that this coherence is produced by the radiative coupling between two 
atoms even when the dipole matrix elements d\ and d,2 are orthogonal. 
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FIG. 5. Time evolution of the coherence in the excited states of atom A with atom B being in |3s). The coherence evolves 
only for non-zero values of r„ c and ft vc . Large oscillations are seen in p\ 2 which is decided by Q vc , and S. The parameters 
used are 9 — tt/2, (f> = 7r/4, R = \/2ir and 8 — 3j. 

In Figs. [6-7], we plot the probabilities that atom A is in state \%a) and atom B is in \jb), which we denote by 
Pi-j = (iA\(jB\p(t)\iA)\jB}- In Fig- 6(a), we present P3;2(i) that represents the simultaneous probability of atom A 
being de-excited to state \3a) and atom B being excited to the state \2b)- Fig. 6(b) is the plot of P2;3(t) that represents 
the probability that the atom A is excited to state \2a) with atom B being in |3b). Obviously both p^ and P2-3 
become zero if T vc = D, vc = 0. It is observed that smaller the atomic separation larger is the excitation probability. 
For atomic separation R = X/2n, the excitation probabilities are very large, e.g. more than 25% of the population in 
atom B could be excited to state \2g) at t ~ O.3/7 (Fig. 6(a)) and, similarly in atom A, ~ 18.5% population could be 
excited to state \2a) at t ~ O.5/7. Thus significant amount of energy transfer can take place between the states 1 1^4.) 
and \2b), though the corresponding transition dipoles are orthogonal to each other. Note that the initial evolution of 
P2;3 is much slower compared to the evolution of p^p. This can be understood as follows: The excitation of atom B 
to the state \2b) can be caused by a single photon transfer from A to B [the process |1a,3b) — > |3a,2_b)], whereas 
the excitation of atom A to the state 12^) occurs only through atom B and this involves a net transfer of two photons 
[processes |1a,3_b) — > 13,4,2b) — > |2a, 3s) or |l^,3s) — > \3aAb) — > |2a,3_b)]- The oscillatory character of p^-2 and 
P2;3 comes from non-vanishing 5 and from the dd-coupling co-efficients £li and £l vc . The excitation probabilities are 
seen to be larger for degenerate excited states (S — 0) compared to that with finite separation between the excited 
states. For very large S (>> 7), this interference effect disappears. 
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FIG. 6. The time evolution of (a) P3-2 and (b) p2 ; 3 are plotted for different parameters with the initial condition that atom 
A is in state |1a) and atom B in |3_g). In both cases 9 = tt/2 and <j) = 7r/4. The values of other parameters are shown as 
legends. 



In Fig. 7, we present a comparative study of the probability that atoms remain in their initial states, i.e. pi-a, in 
the presence and absence of the dd-coupling terms. The probability of atom A staying in decays exponentially in 
the absence of atom B. However, in the presence of the second atom, the nature of its decay is significantly modified - 
large oscillations are seen in p 1; 3 in the presence of the new coherence terms, which is evident from Fig. 7. The origin 
of this oscillation is attributed to the large values of fi„ c . 
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FIG. 7. Plot of the probability that both atoms remain in their initial states. This probability is plotted on a logarithmic 
scale as a function of time in linear scale. The parameters used are 9 — it/2, (f> = 7r/4, R = A/4 and 5 = 0. 



V. TWO ^-SYSTEMS WITH MAGNETIC SUB-LEVELS IN THE PRESENCE OF A MAGNETIC FIELD 

In this section we consider the new coherence effects in two U-systems with m-degenerate magnetic sub-levels as 
excited states. The system could be, for example, a 40 Ca system - where 4 1 P\ degenerate sub-levels would correspond 
to the excited states = \ j = l,m = 1) and |2 M ) = \j = l,m = —1), and the 4 1 S'o state would correspond to the 
ground state |3 M ). In this case the dipole matrix elements d\ and d 2 are complex and orthogonal to each other, and 
are given by 

-* X zb zy 

di = -de-, d 2 = di + ; e± = — ^=-, (28) 

V2 

where d is the reduced dipole matrix element. The magnetic field produces a Zeeman splitting 5 and fixes the 
quantization axis (z axis in our case). The geometry can be taken to be the same as in Fig. 2. However, in the present 
case, d\ and d 2 being complex dipoles, they are not fixed along the real axes unlike in Fig. 2. Using Eq. (^l|), the dd 
coupling coefficients for this scheme can be obtained 

Tl = T 2 = M-Im ( X xx + Xyy) = ^( 2P * - Shl2 
n i= n 2 = i|-Re {Xxx + Xyy) = ^( 2P r ~ ^ 9Q r ), 

r KC = [Im (Xxx - Xyy) + i Im (Xx y + Xyx)} = ^ sin 2 9e 2 ^Q l , 

V vc = [Re (xxx - Xyy) + i Re ( X xy + Xyx)] = ^ sin 2 9e 2 ^Q r . (29) 

The Ps' and Qs' are as defined in Eq. (26). In deriving (p9|), we have used the fact that Xxy = Xyx- It may be 
noted that Ti and £li are real, and are independent of the azimuthal angle, whereas r„ c and Q vc are complex and 
are functions of 4>- Fot 9 = nir, the coherence terms disappear in Eq. (|i~8|). Thus if R is perpendicular to the plane 
containing both the dipoles, i.e. both atoms lie on the quantization axis (z-axis), the coherence effects vanish. 

The solutions of the master equation can be recalculated using the above coefficients and the analog of all the 
results presented in Sec. IV can be produced for the present system. For completeness, we present the numerical plot 
that shows the excitation probability ps ]2 with the initial condition pis{t = 0) = 1. The time evolution of p% [2 is 
similar to the one in the case of real dipoles (cf. Fig. 6(a)). 
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FIG. 8. The time evolution of the probability p-^-2 when the dipoles di and d,2 are complex. All the parameters are same as 
in Fig. 6. 

It may further be noted that P3 ; 2 is independent of <fi though T vc and £l vc are functions of cf>. This is because P3 ; 2 
is a function of the absolute values of r. uc and fl vc . - which can be shown from ( |25| ) and (^9|), to the lowest order in 
r„ c and f2 uc , 




We now discuss how the new coherence effect can be monitored experimentally for the above mentioned system. 
The dipole transitions <-> |3 M ), in the system described above, involve photons having <r + polarization. Thus the 
emission from 1 1 M ) — ► |3 M ) does not contain any field component in er_ polarization. On the other hand, the emission 
from |2 M ) — > \3fj) would contain cr_ component. Thus the signal that one has to look for is - the intensity of the 
emitted photon from |2 M ) levels in er_ polarization, which would a be measure of the total excitation probability to 
|2 M ) states and hence would confirm the occurrence of VIC. Another possibility to probe the population in |2 M ) will 
be to excite it with a circularly polarized radiation to a fourth state 6 1 Sq and to monitor the fluorescence from 6 1 Sq. 

VI. CONCLUSIONS 

In conclusion, we have shown that the radiative coupling between multilevel atoms with near-degenerate transitions 
can produce new interference effects which are especially important when the distance between two dipoles is less 
than a wavelength. We have demonstrated this possibility by considering two identical V-systems such that the pair 
of transition dipole matrix elements in each system are orthogonal to each other in both the atoms. Such interference 
effects are especially significant in the energy transfer studies. The choice of orthogonal dipole matrix elements enables 
us to specially isolate the effects of the vacuum induced coherences in the radiative coupling between multilevel atoms 
with nearly degenerate transitions. We have presented detailed numerical results to bring out the role of multi-atom 
multilevel interference effects. 
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